Comment for physicists, written 30 January 2011 

This paper has been published in 1995 in the mathematical journal Advances 
in Applied Clifford Algebras. It was my unique publication on Clifford algebras. 
Nevertheless I was personally invited to participate in several conferences on Clif- 
ford algebras, although I am a physicist, and a use of Clifford algebras is only a 
mathematical method of the Dirac equation representation. In my paper the 7- 
matrices are eliminated, and the Dirac equation is presented in the hydrodynamical 
form, where the 7-matrices are absent. On one hand, the physicists do not use the 
Sauter - Sommerfeld method of the 7-matrices representation in the form of Clif- 
ford algebras. On the other hand, one can eliminate 7-matrices, only if they are 
not represented in a concrete representation (i.e. they are to be considered as basic 
elements of a Clifford algebra). 

Elimination of 7-matrices admits one to obtain essential physical results: (1) 
discovery that the Dirac particle has a complicate internal structure (it is not a 
pointlike particle), (2) discovery that the internal degrees of freedom are described 
non-relativistically. These statements are very essential, but they are rather un- 
expected simultaneously, and one cannot understand them without a knowledge of 
the Sauter - Sommerfeld method. I tried to publish the physically essential parts 
of the paper [U [2] in physical reviewed journals. Unfortunately, I failed, because, 
apparently, the reviewers did not know the Sauter - Sommerfeld method of the 
7-matrices representation. They do not trust in mathematical correctness of my re- 
sults, although experts in Clifford algebras, publishing me paper in the mathematical 
journal did. 
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Abstract 

The distributed system Sd described by the Dirac equation is investigated 
simply as a dynamic system, i.e. without usage of quantum principles. The 
Dirac equation is described in terms of hydrodynamic variables: 4- flux j l , 
pseudo- vector of the spin S l , an action tup and a pseudo-scalar k. In the 
quasi-uniform approximation, when all transversal derivatives (orthogonal to 
the flux vector j l ) are small, the system Sd turns to a statistical ensemble 
of classical concentrated systems Sdc- Under some conditions the classical 
system Sdc describes a classical pointlike particle moving in a given electro- 
magnetic field. In general, the world line of the particle is a helix, even if the 
electromagnetic field is absent. Both dynamic systems Sd and Sdc appear to 
be non-relativistic in the sense that the dynamic equations written in terms 
of hydrodynamic variables are not relativistically covariant with respect to 
them, although all dynamic variables are tensors or pseudo-tensors. They be- 
comes relativistically covariant only after addition of a constant unit timelike 
vector p which should be considered as a dynamic variable describing a space- 
time property. This "constant" variable arises instead of 7-matrices which are 
removed by means of zero divizors in the course of the transformation to hy- 
drodynamic variables. It is possible to separate out dynamic variables k, k 1 
responsible for quantum effects. It means that, setting k, k 1 = 0, the dynamic 
system Sd described by the Dirac equation turns to a statistical ensemble 
£Dqu of classical dynamic systems Sdc- 

1 Introduction 

The Dirac equation considered as a dynamic equation for a wave function in frames 
of quantum mechanics principles has been investigated almost completely, and it is 
hardly possible to add anything new. At the same time the Dirac equation consid- 
ered as a dynamic equation for a distributed dynamic system Sd displays a series 
of such unexpected properties as existence of dynamic variables k responsible for 
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quantum effects and appearance of a constant timelike vector describing a split of 
the space-time into the space and the time. 

From mathematical point of view our investigation of the Dirac equation is simply 
a change of variables, when the Dirac four-component complex wave function is sub- 
stituted by certain tensor variables having also eight independent real components. 
The tensor variables are: the 4-flux vector j l , i = 0, 1,2,3, the spin pseudo- vector 
S l , i — 0,1, 2, 3, the action scalar <p, and a pseudo-scalar k. These quantities will be 
referred to as hydrodynamic variables. 

From physical standpoint the dynamic system Sd is considered as a distributed 
dynamic system describing in some way an electron (or positron) motion. The mo- 
tion of a single particle is supposed to be stochastic, the variables j % describing the 
mean 4-flux of particles. Interpretation of other dynamic variables of Sd is produced 
on the base of a comparison with a statistical ensemble of classical particles. The 
quantal correspondence principle, when a linear operator corresponds to any phys- 
ical quantity, is not used. Other quantum principles are not used also. Instead of 
them one uses a more general statistical principle which asserts: A set S (statistical 
ensemble) of many similar independent stochastic systems S s is a deterministic dy- 
namic system Sd [1]. The term "a stochastic system" means that experiments with 
a single system are irreproducible, and there are no dynamic equations for such a 
system. If S means either stochastic system S s , or a deterministic system Sd, the 
statistical principle can be formulated in the form 

S[S] is S d 

The statistical principle can be applied both to quantum and classical systems. It is 
a more general statement, than the set of quantum principles. In particular, it can 
be valid for such stochastic systems which are neither quantum, nor classical. The 
main concept of this approach is the statistical ensemble considered as a dynamic 
system (not a probability density, or a probability amplitude). This approach will 
be referred to as a statistical ensemble technique (SET). 

There is a reason for consideration of the Dirac equation without usage of quan- 
tum principles, i.e. simply as a dynamic equation for the system Sd- Recently one 
suggested a hypothesis that the real space-time is a distorted Minkowski space-time, 
and the distortion of the space-time is a reason for quantum effects [2]. Distortion is 
such a deformation of the Minkowski space-time which transform one-dimensional 
world lines into three-dimensional world tubes. The world tubes of particles appear 
stochastic, and their statistical description coincides with quantum description pro- 
vided the Planck constant h determines the space-time distortion (thickness of the 
tubes). This hypothesis was proved for a non-relativistic free particle. It is interest- 
ing to test the hypothesis in the case of a relativistic particle. For such a test it is 
necessary to separate the action describing the ensemble of Dirac particles into two 
parts: a classical part and a part responsible for quantum effects. 

Another motive is as follows. The conventional quantum dynamics technique 
(QDT) is an axiomatic construction. QDT relates to SET approximately in the 
same way, as the axiomatic thermodynamics relates to the statistical physics. For 
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instance, the Brownian motion cannot be explained and understood from standpoint 
of thermodynamics. Thermodynamic axioms do not permit to do this. The statisti- 
cal physics explains the Brownian motion phenomenon and restricts an application 
of thermodynamic principles, because a statistical approach is more general, than 
the approach of the axiomatic thermodynamics. 

Something like that one can see in the field of the quantum dynamics. Problem 
of pair production is a principal problem of the high energy physics. In those areas, 
where the pair production is unessential and can be considered as a correction (for 
instance, in quantum electrodynamics) the quantum theory succeeds. In those areas, 
where pair production is a dominating effect, the quantum theory (in particular, 
QFT) failed. 

For last fifty years the quantum field theory has not succeeded in solving the 
problem of pair production. A suspicion arises that this failure is not accidental. 
Maybe, the problem of pair production cannot be solved in framework of the quan- 
tum principles, and a more general statistical approach is necessary. 

Classical limit of the Dirac equation is rather difficult to obtain, because it 
contains such non-classical quantities as Dirac 7-matrices, which hardly can be con- 
sidered from the classical point of view. 

Usually a quantum electron moving in a given electromagnetic field Ai, I — 
0,1,2,3 is described by the Dirac equation 

-ih^di^ + eAi^ijj + mijj^O, (1.1) 

where ip is a four-component complex wave function. The light speed c is chosen to 
be equal to 1. It is possible to transform the variables ip and to describe this system 
in terms of variables if, j l , S l , (/ = 0, 1, 2, 3), k, defined by the relations 

y = ^ 7 V, / = o,i,2,3, v> = vV; 

S l = ^7 5 7V, J = 0, 1, 2, 3, 75 = 7° 123 = 7V7V; 

dup = - W)(2n^)-\ 1 = 0,1, 2, 3; (1.2) 

cos/t = ^(j l ji)~ 1/2 ; 

Here 7', / = 0, 1, 2, 3 are the Dirac 7-matrices satisfying the commutation relation 

7*7* + 7*7* = 20* i,k — 0, 1, 2, 3, (1.3) 

where g lk =diag(l, — 1, —1, —1) is the metric tensor. Only two components of the 
pseudo- vector S l are independent, because there are two identities 

S l Si = -j l j h j l S l = 0. (1.4) 

Pseudo- vector S l is treated as a spin pseudo-vector, because it is connected uniquely 
with the spin tensor S ml ' h , defined by the relation [3] 

S mi,k = l^( 7 V m + <r im 7 fc )^, a lm = ^(7*7 m - 7"V) (1.5) 



4 



The relation between S l and S ml ' k has the form 

S i = ±g ii e jmlk S ml > k , i = 0,1,2,3; 

S ml > k = ~e imlk S h m,l,k — 0, 1, 2, 3. (1.6) 

Here Sj m ik and e %mlk are Levi-Chivita pseudo-tensors (^0123 — 1, e 0123 = —1). 

Thus, the scalar </?, pseudo-scalar /t, vector j', and pseudo- vector S l have 8 real 
independent components which are used instead of 8 real components (4 complex 
components) of the Dirac wave function if). 

Realization of such a transformation is a goal of the present paper. Such a 
description can be regarded as a description in terms of hydrodynamic variables. 
Under some conditions the dynamic system Sp turns into a statistical ensemble 
&Dqu of classical dynamic systems Sdc- The Sdc can be interpreted as a classical 
analog of the Dirac electron. £oqu is described usually in terms of hydrodynamic 
variables j l , ip, The dynamic systems Sp and Soqu are distinguished by some 
terms of their Lagrangians, but not conceptually. Thus, one can consider a dynamic 
distinction between the Dirac dynamic system Sjj and the corresponding classical 
ensemble. 

In the second section some properties of the statistical ensemble E c \ are investi- 
gated. The third section is devoted to introduction of hydrodynamic variables. In 
the fourth section the action for the Dirac equation is written in terms of hydrody- 
namic variables. The fifth section is devoted to consideration of quasi-uniform state 
of the system Sd- Relativistic invariance of dynamic equation in terms of hydrody- 
namic variables is discussed in the sixth section. The seventh section is devoted to 
investigation of the classical analog Sdc of the Dirac electron. 



2 Statistical ensemble of classical dynamic 
systems 

Let there be a classical system S described by the Lagrangian function L(t, x, dx/ dt), 
where x = {x a }, dx/dt = {dx a /dt}, a — 1,2, . . .n are generalized coordinates and 
velocities. Then, by definition a pure statistical ensemble £ c i of systems S is a set 
of similar independent systems S. Its Lagrangian is a sum (integral) of Lagrangians 
L. The action for the ensemble £ c i has the form 

/ct=n 
L(t, x, dx/dt)dtd£, d£ = JJ d£ a (2.1) 

a=l 

where £ = a = 1, 2, ... n are Lagrangian coordinates labelling the systems of 

the ensemble, and x = x(t, £) is a function of t and 

The dynamic system £ c i can be considered as a fluid described in the Lagrangian 
coordinates (the time t and Lagrangian coordinates £ are independent variables). 
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The same action written in the Euler coordinates (t, x are independent variables) 
has the form 



/n 
{C(x,j) - hf[d lV + 9 a (mU}d n+1 x, d n+1 x = J] dx\ (2.2) 

where j = {j°, j}, j = {j a }, a = 1, 2, . . . n, tp, £ = {£ Q }, a = 1, 2, ... n are functions 
of x = {a; ,x}, x° = t, x = {rc Q }, a = l,2,...n. A summation is made over 
repeating indices: over Latin ones (0 — n) and over Greek ones (1 — n). All variables 
j = v 9 ; £ are functions of x. The function C(x,j) is defined by the relation 

C(xJ)=?L{aP,x,i/j°) (2.3) 

where L is the Lagrangian of a single system S. The variables j = {j°, j} describe 
a flux of particles in the (n + l)-dimensional space V of coordinates x. They are 
connected with the variables x(t, £), dx.(t,£)/dt by means of relations 

J'° = det||^||, ^ = 9^, a,/3 = l,2,...n, j=J°^ (2-4) 

and the functions £ = ^(x) are determined as solutions of the equations 

x a = x a (t,$,), a = l,2,...n. 

Functions g a (£) are arbitrary functions of Lagrangian coordinates £ a , a = 1, 2, ... n. 
Appearance of arbitrary functions is a result of integration of some dynamic equa- 
tions arising in the course of transformation from the action (2.1) to the action 
(2.2) (see Appendix A). In turn a possibility of such an integration is connected 
with an invariance of dynamic equation with respect to arbitrary transformation of 
Lagrangian coordinates £ 

= !«(€), det || dlp/dU ||= 1, a, = 1,2,... (2.5) 

Flux j = {j°, j} is invariant with respect to the renumbering transformation (2.5). 

A statistical ensemble of classical systems is described conventionally by the 
distribution function F(t, x, p) in the phase space of coordinates x and momenta p. 
For the pure statistical ensemble the distribution function has a special form 

F(t,x,p)=/(t,x)5[x-P(t,x)] (2.6) 

where 

P = {Pa}, Pa = dL/d(dx a /dt), a = l,2,...n (2.7) 

P(t, x) = {P a (t, x)}, a — 1, 2, . . . n is a set of functions depending only on t and x, 
and 5 denotes the Dirac 5-function. Evolution of this ensemble is described by the 
Liouville equation of the form 

dF dH dF dH dF _ 
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where H = H(t, x, p) is the Hamiltonian function of the dynamic system S. Dy- 
namic equations for the variables j°, P(t,x) can be obtained by means of a substi- 
tution of Eq. (2.6) into Eq.(2.8). 

Among three ways [(2.1), (2. 2) and (2.8)] describing the pure statistical ensemble, 
the action (2.2) is most convenient for comparison with action of the system Sd- 

Dynamic equations, determined by the action (2.2), have the form 



^ = m? = o, 

dip 



5A 



dL 

— - mp - 



o, 



0,1, ...n; 



6A 



-h 



d^a die 



3%U = 0, 



a 



.n; 



(2.9) 

(2.10) 
(2.11) 



j l , I — 0, 1, . . . n is the current (n + l)-vector in the space V. The vector field f is 
tangent to the trajectories of systems in V. According to Eq.(2.3) the vector 



Pi = 0ji' Z = 0,l,...ra; 



L — 



dL 



d(dx a /dt) 



dx a /dt=j a /j° 



Pa 



dL 



d{dx a /dt) 



a = l,2,...n (2.12) 

dx a /dt=j a /j° 

associates with the canonical momentum of a single system of the ensemble. Thus, 
according to Eq.(2.10) 



Pi = h(d i <p + g a (S)t a J, i = 0,1,2,... n 



(2.13) 



where h is a constant having a dimensionality of the action. In this case tp, g a , and 
^ can be considered as dimensionless quantities. H can be treated as the Planck 
constant, although in the given case it has no quantum meaning. 

Eliminating variables ip and £ from Eqs. (2.10), (2.11), one obtains the equations 



mp k - dm\ = - 9 k d -%^\ - o. 



Of 



dji 



k = 0,l,...n (2.14) 



Together with Eq.(2.9) the equations (2.14) form a system of n+1 dynamic equations 
for some kind of a fluid described by the current vector j l . 

At the linear transformation of coordinates x l , i — 0, 1, . . . n the current vector 
j l , I = 0, 1, . . .n transforms as a vector. In this case the Lagrangian coordinates 
transform as scalars. But the consideration of £ as scalars is rather conventional, 
because of the renumbering transformation (2.5). 

At the transformation (2.5) one obtains 



(2.15) 
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In other words, g a (£) transform at the renumbering transformation as components 
of a vector in the space of Combining any linear transformation of coordinates x l 
with some renumbering transformation, one can ascribe practically arbitrary trans- 
formation properties to the Lagrangian coordinates 

Any renumbering transformation (2.5) is a kind of a gauge transformation, be- 
cause the renumbering changes a description of the state of the statistical ensemble 
without changing the state itself. It is easy to verify that a set of the renumbering 
transformations forms a group. 

For the pure statistical ensemble of classical pointlike charged particles moving 
in the given electromagnetic field the action (2.2) takes the form 

Mj^M = J [-rnVMi-eA i j l -h?{a i <p + g a (t)Z 0h J]d i x (2.16) 

where the speed of the light c = 1, and summation is made over repeated indices (0 
- 3) for the Latin indices and (1 - 3) for the Greek ones. 



3 Transformation of variables 

Transforming the Dirac equation (1.1) to the new variables (1.2), one uses the action 
for the equation (1.1) 



Ad[iP, i>) = j (-mvpvp + -hi^ydiijj - -hdiipj ip - eAiipj ip)d x (3.1) 

Expressing the variables (1.2) through the wave function if>, one uses the technique, 
where the wave function is considered as a Clifford number with 16 base units: 
I, 7 l , 7 lfc , 7 lfc ', j lklm (all indices are different, and Clifford numbers satisfy (jl.3jl ). 
Reduction of the Clifford numbers is realized by means of zero divisors [4,5]. 
Let us introduce matrices 7 5 , cr = {a a }, ot = 1, 2, 3 

7 5 = 7° 123 = tVtV, <xi = -*7 23 , a 2 = -* 7 3 \ « 3 = -*l 12 , 

l 5 (?a = Va7 5 , 7°° = -ils^w a = 1,2, 3; (3.2) 
7 °cr = <T7°, 7°7 5 = -7 5 7° 
According to Eqs. (1.3), (3.2) the matrices cr = {a a }, a = 1, 2, 3 satisfy the relation 

Cq "/? = $0/3 + ^a^7 cr 7, a, /3 = 1,2,3 (3.3) 

where is the antisymmetric pseudo-tensor of Levi-Chivita {eva = 1)- 
Let us define the wave function tp in the form 

i> = n^*7°, ip* = AUe-^e-^'^e-^^ (3.4) 
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where (*) means the Hermitian conjugation, and 



n = i(l + 7°)(l + «r) 



z = {z a } = const, 



a = 1,2, 3; 



1 (3.5) 



is a zero divisor. The quantities A, k, <p, r] — {r] a }, £ = {C*}, & = 1,2,3 are nine 
real parameters, determining the wave function ip. 

Using relations (3.2), (3.3), (3.5), it is easy to verify that 



it = n, 



7 n = n, 
na„n = z a u. 



zo-n = n, n 7s n = o, 

1,2,3. 



a 



(3.6) 



Generally, the wave functions defined by Eq.(3.4) are 4x4 complex matrices. 
In the proper representation, where II has the form 



IT 



/ 1 








o\ 


























1° 








0/ 



(3.7) 



the ijj* have the form 



Their product ip*Oip has the form 



(^1 








o\ 













^3 











V ^4 











/ r 2 r 3 i>\ \ 





\ / 



(3. 



( a \ 





\ J 



all = Ua 



(3.9) 



where O is an arbitrary 4x4 matrix and a is a complex quantity. If / is an analytical 
function having the property f(0) = 0, then the function /(all) of a 4 x 4 matrix 
of the type (3.9) is a matrix /(a)II of the same type. For this reason one will not 
distinguish between the complex quantity a and the complex 4x4 matrix all. In 
the final expressions of the type all (a is a complex quantity) the multiplier II will 
be omitted. 

By means of relations (3.2), (3.6) one can reduce any Clifford number IlOn to 
the form (3.9), without using any concrete form of the 7-matrix representation. This 
property will be used in our calculations. 

By means of H3.4I) the variables ?/>?/>, j , S l , I = 0, 1,2,3 defined by expressions 
(1.2) can be presented in the form 



W = ip*j°ijj = A 2 Ile 75K n = A 2 cos fill 



(3.10) 
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j°n = VtV = A^e-^e^'^e^Il = A^e-^m = A 2 cosh^n (3.11) 
where 

S = {S 1 ,S 2 ,S 3 }, E a = e-3^(T Q e^ a = 1,2, 3; (3.12) 

The matrix S Q satisfy the same relations (3.3), as a a do. 
In the same way one obtains 

j a u = ^* 7 0a vn = A 2 ne~t^ s?7 (-i7 5 S a )e-t 7 5 ST 'n = 

= A 2 n(coshf -i 75 Svsinhf)(-i 75 S a )(coshf -i 75 Svsinh§)n = (3.13) 
= A 2 smh(ri)v a Il, a = 1,2,3, 



where 



v = {</*}, v a =r] a /r], a = 1,2,3; v 2 = 1. (3.14) 



s°n = ?/>*H7 5 )v = A 2 n(-i 75 ) e -^s»7n = 

= A 2 U sinh(??)£vIT = A 2 sinh^IT, 

S a U = ip*a a if;U = A^e^^ei^m = 

A 2 [C + (cosh??- l)v a (v£)]Il, a = 1,2,3. 

Here £ = a = 1, 2, 3, are determined by the relation 



(3.15) 



f*n = IT£ a n, a = 1,2,3 (3.16) 
It follows from Eqs.(3.11), (3.13) 

f 3l U = A 4 U, A={j l 3l ) 1 l^p 1 / 2 (3.17) 
According to Eqs.(3.12), (3.16) one obtains 

£11 = {[z - n(nz)]cosC + (z x n) sinC + n(nz)}II (3.18) 

where 



£ = {?}, a = 1,2,3; (^C 2 ^, 

n = C/C, n 2 = l (3.19) 

The wave function (3.4) depends on 9 real parameters: A, if, k, i] a , ( a , a = 1,2,3. 
The wave function has 8 real independent components, and not all parameters 
A, <f, k, i] a , ( a , a — 1, 2, 3 are independent. Let us fix one of the parameters, namely 
let us set 

C = (3.20) 

Then Eq.(3.18) takes the form 

£ = 2n(nz)-z (3.21) 
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It can be solved with respect to n = One obtains 

C/vr = n =(£ + z)[2(l + zOr 1/2 (3.22) 

Using Eqs.(3.11), (3.13), (3.17), (3.22), one can express parameters A,ip, K,i] a ,( a , 
a = 1,2,3 describing the wave function, through the variables j l , S l ,ip, k. One 
obtains 

-lroa _ J a ( S 3 ) i n = 10.Z- n=,fJ7, 



C = f ^[S a - f {f l p) ^ « = 1.2,3; p=y/H (3.23) 
cosh V =j°/p, v a = n — = -= £ , , a = 1,2,3 (3.24) 

v v ) -j ji 

Using Eqs.(3.4), (3.20) - (3.24), one can present the wave function (3.4) in terms of 
variables j\ S\ k, ip: 

* = 2V /(1 + 0(J° + P) :[(J ° + + ^ " ll5S0 ~ ' 75(j + ^ X S)CT]I1 (3 ' 25) 

where £ is expressed through j l ,S l by means of Eq.(3.23), the symbol x means a 
vector product, and II is defined by Eq.(3.5). 

Any expression of the form ipOip, where O is an arbitrary combination of 7- 
matrices, can be expressed through the variables S\j\ k. For instance, 

^(y 7 * _ 7 *y^ = y p ^" sl - ^) sin K + £lkim ji s m cos «]n. (3.26) 
4 Transformation of the action 

Now let us calculate the expression (two middle terms of Eq.(3.1)) 

+ h.c. = \hi? [(d - i 75 <rV) (i(p + i 75 /c)] iJ + h.c. 

+iM 2 ne~5^ e -^75^( ( 9 _ i 75< rV)(e-^<"7 e §^)n + h.c. 

(4.1) 

where "h.c." means the term obtained from the previous one by the Hermitian 
conjugation. Using relations (3.11)-(3.15), the expression (4.1) reduces to the form 

|fa/>7 Z <9;-?/> + h.c. = —hj l di(f — ^hS l diK 

+ fM 2 IIe-^ ST '(<9o - i7 5 £V)e^5 Sj ?n + h.c. 
+fM 2 ne- i7 5 s ^e-^<9 ei' TC II + h.c. 

+ | M 2 ne"t^Sr,(_ i75Sa ) e -f75Sr ?e -f<TC 9Qe f ( TCn + h.C. 

(4.2) 

where the matrix X is not differentiated. 

Taking into account the fourth Eq.(3.6) and Eqs. (3.11), (3.13), the expression 
(4.2) reduces to the form 

\h^ l dii) + h.c. = -hj l di(p - \hS l diK + ^hj l Ue-^die^ c U + h.c. 

+ fM 2 ne-^ 5 ^(9 -i7 5 SV)e-^ s '?n + h.c. { ' } 
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Substituting the relation (3.20) into the third term of Eq.(4.3), one obtains by means 
of Eq.(3.22) 

^-HfUe-^die^U + h.c. = ^-hj l Ua a a ^riad^U + h.c. 



2(1 +£z 

where 

g a {$,) = - l -e aM ez < {l+^r 1 (4.5) 
Calculation of the last term of Eq.(4.3) leads to the following result 

F 4 = {hA 2 He-^^(d Q - ry 5 EV)e-^5ST7n + h.c. = , . 

-\hA 2 e aM [d a rivPC + sinhr^^ 7 + 2 sinh 2 (2)^d A 7 ]n 1 j 

Let us introduce two constant vectors 

f = {1,0,0,0,}, z l = {0,z\z 2 ,z 3 } (4.7) 

which satisfy the following conditions 

/7j = 1, A = 0, A = -l. (4.8) 

By means of relations (3.11), (3.13) and (4.7) the expression (4.6) reduces to the 
form 

= - 2{p ^ fSjs f ki m [d k (f + fp)](j l + fp)[C - f m (Cfs)] (4.9) 

where Sikim is the Levi-Chivita pseudo-tensor (£0123 = 1) and £ m = The 
value of £° is unessential. 

Let us introduce the unit timelike vector 



ji + f i p ji + fi p 



(4.10) 



vu i +f i p)Ui+fip) VMp+m 

and two spacelike vectors 

/''' C- f'lCl). < = 0,1,2,3; ^ = -1, (4.11) 

H = , = . = . 4.12 

Then according to Eqs.(4.3), (4.4), (4.9) the relation (4.2) can be presented in the 
form 

'-h^drt + h.c. = -hj l [d^ + 9a(mC] - \s l d t K + hpe tklm q\d k q l )v m (4.13) 
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where the second term can be written also in the covariant form 

- frj l 9a(me = -frj l e jkS mPi j (dipi k )f s z m (4.14) 

Now by means of relations (3.10)-(3.12) one can present the action (3.1) in the 
form 

A D [j, <p, k, f] = J (£ d + C ql + C q2 )d 4 x (4.15) 

C cl = -mp - eAa 1 - hf[d lV + 9a(mC], P = vWi (4.16) 
£ ql = 2mpsm 2 (^)-^S l d l K, (4.17) 

C q2 = hpe lklm q\d k q l )u m (4.18) 

g a (£) is defined by the equation (4.5). S l , I = 0, 1,2,3 are considered as functions 
of j l and f , defined by the relations 

S°=te, S a = pC + ^ I i T , a = 1,2,3 (4.19) 

p + rfk 

obtained from Eqs (1.4), (3.23). Not all variables £ = ot = 1, 2, 3 are indepen- 

dent, because they satisfy the restriction 

£ 2 = TC = 1 (4-20) 

as it follows from Eqs. (3.5), (3.19), (3.21). Variation of the action (4.15) with respect 
to £ a , a = 1, 2, 3 under the condition (4.20) leads to the dynamic equations 

^(^-^ = 0, = 1,2,3 (4.21) 

There are only two independent equations among the equations (4.21), because a 
contraction of Eq.(4.21) with leads to an identity. Note that the term (4.16) of 
the action (4.15) coincides with the Lagrangian of the action (2.16) for the statistical 
ensemble of classical pointlike particles. 

Eliminating k from the action (4.15) and the dynamic equation SAd/Sk, = 0, 
one can write the action (4.15) in the form 

An\j, <P,t]=J (4* + CT ql + C q2 )d A x (4.22) 

/ ft 2 ft ftr) <?' 

4* + C? ql = -yjmWi - ^(d^r - fa arcsin (^)- 

- ej% ~ hf[d tV + g a (mC} (4-23) 

S\ i — 0, 1,2,3 are functions of £ and j l which are determined by Eqs. (4. 19). C q2 
is determined by Eq. (4.18). 
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5 Classical Part of the Action 



The classical part of the action Ad can be separated out either by vanishing h, or 
by a usage of the quasi- uniform state of the system, when all spatial gradients are 
small, and quantum effects disappear. 

In the non-relativistic case a quasi-uniform state satisfies the condition 

I— Vu|«|u|, u = f,j 1 ,f,f,K. (5.1) 
m 

which is written in the coordinate system, where | j \<g\ j° |. In the general case 

such a coordinate system does not exist, and the condition (5.1) is written in the 
form 

h 

\ —rd iU \<^\ u \, u = j ,j\j 2 ,f,K, (5.2) 
m 

where /* is any unit vector orthogonal to j l 

1% = 0, 1% = -1. (5.3) 

It means that all derivatives across the direction of the vector j l are small. 
Let us represent all derivatives in the Eq.(4.18) in the form 

d k q l = d k ± q l + J -^-d a q l , d k ± = d k - 3 -^-d s (5.4) 
p z p z 

The first term in rhs of the first equation (5.4) describes a transversal part of the 
derivative, i.e. a derivative in the direction orthogonal to the vector f . 

According to Eq.(5.2) all terms containing the transversal derivative d\ are small 
with respect to the first term of Eq.(4.16). Indeed, an estimation of the transversal 
part of the Lagrangian (4.18) has the form 

3 

I Ti P E; Mm q\d\q l )v m |^ hp | J2 l U 9 ^ l« m P> ( 5 - 5 ) 

s=0 

where ZK, s = 0, 1, 2, 3 are unit vectors orthogonal to f 



l l - jsrJ^iZP^Ml s - n 1 2 3 (5 6) 

(no summation over s) 
Neglecting the transversal part and using Eqs.(4.10), (4.11), one can write Eq.(4.18) 
in the form 

C qu2 = hp- x 2 s E lklm q l 3\d s q l )v m = . ^ == e iklm fj k d s q l C 
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The relation (4.17) can be written in the form 

k h ■ ■ S i 

C a i = mp(2sin 2 w l &re), w l = , w*7j = 0. (5.8) 

9 2 2m y ' y/-S l Si 

Let us take into account that the last term in the first equation (5.8) is small 
with respect to the first term of Eq.(4.16). Then neglecting small terms and taking 
into account Eq.(5.7), one obtains for the action (4.15)-(4.18) in the quasi-uniform 

state 

A Dqu [j, ip, re, £] = J {-mp cos re - eAif - hf\d^ + g a (£)di£ a ] 

where g a {£) is determined by Eq.(4.5). 

Total derivatives p~ l j l di of variables j l , re, £ in the dynamic system Soqu de- 
scribed by the action (5.9) are determined mainly by dynamic equations. Derivatives 
in the orthogonal directions are determined by the initial conditions which must be 
such ones, that these derivatives were small enough. 

Although the action (5.9) contains a non-classical variable re, in fact this variable 
is a constant. Indeed, a variation with respect to re leads to the dynamic equation 

SAoqu _ m ^ sin re = (5.10) 
ore 

which has solutions 

re = nn, n = integer (5-H) 
Thus, the effective mass m e ff = m cos re has two values 

m e ff = m cos re = ±m (5-12) 

The value m e ff = m > 0, (re = |mr) corresponds to a minimum of the action 
(5.9), whereas the value m e ff = — m < corresponds to a maximum. Apparently, 
m e ff > corresponds to a stable ensemble state, and m e ff < does to unstable 
state. 

Eliminating re by means of the substitution k = ^nir in Eq.(5.9), one obtains the 
action 

A Dqu [j, (p,£] = J {-mp - eAif - hf[diip + g a (€)diC] 

+ TT^TTeHwWr}^ (5-13) 

Let us introduce Lagrangian coordinates r = {t^}, i = 0,1,2,3 by means of 
relations 

.. 8D _ d{x\T 1 ,T 2 ,T 3 ) _ d(T ,T 1 ,T 2 ,T 3 ) _ 



5ro,i (9(x°, x 1 , x 2 , re 3 ) <9(x°, x 1 , x 2 , x 3 ) 
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T k ,i = diT k , i,k — 0, 1,2,3 

Taking into account that 

D j diU = D diU = — r = — (5.15) 

OTo,i o{to,t 1 ,t 2 ,t 3 ) dr 

d A x = D~ l d 4 T = D- x dr Q dT (5.16) 

o(x u ,x L ,x 2 ,x iS ) 

the action (5.13) can be rewritten in the Lagrangian coordinates in the form 

AM = j { - m ^-eA, i ' + H^ + H 2V¥ ^ { ^. + . 0) ) 4'r (5.18) 

where the dot means the total derivative d/dr^. x = {x 1 }, i = 0,1,2,3, £ = 
a = 1,2,3 are considered as functions of the Lagrangian coordinates r , 
r = {ri, r 2 , r 3 }. z is a constant unit 3-vector. Ai = Ai(x), i — 0, 1, 2, 3 are function 
of x. The term <9;<p is omitted, because it reduces to a Jacobian (5.17) and does 
not contribute into dynamic equations. 

The action (5.18) describes a statistical ensemble of deterministic dynamic sys- 
tems Sdc- A state of each system Sd c is described by the variables £. The 
variables £ are connected with the spin by the relation (3.23) which takes the form 

(sx) (£x) 
£ = s t== x, s = $ + V ^ ; ^- ~ x (5-19) 

Here s = S/p is the unit spin 3-vector. 

All this means that in the quasi-uniform approximation the dynamic system 
is a statistical ensemble of some deterministic classical systems S dc . The system 
Sdc should be treated as a classical analog of the Dirac electron moving in a given 
electromagnetic field. Under some conditions the classical Dirac electron turns into 
a classical relativistic pointlike particle, but, in general, Sd c is a more complicated 
construction than a pointlike particle. As far as under some conditions Sd is a 
statistical ensemble of classical relativistic particles, one concludes that, generally, 
the Dirac equation describes a statistical ensemble of charged quantum particles (not 
a single particle). Indeed, at first, one knew only that the system relates to an 
electron in some way, but one did not know whether Sd describes a single electron or 
a statistical ensemble of electrons. One discovers that under some (quasi-uniform) 
initial conditions the Sd is a statistical ensemble (of classical systems Sdc)- It means 
that Sd is a statistical ensemble in all other cases. But in the general case Sd cannot 
be an statistical ensemble of deterministic classical systems. It means that Sd is a 
statistical ensemble of stochastic systems. 

As far as the quantum principles are not used, then one uses the statistical 
principle formulated in the sec.l. This principle permits to determine mean values 
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of energy, momentum, angular momentum and other additive quantities for a single 
electron. Indeed, dividing an additive quantity for So by the number of systems 
in the ensemble So, one obtains corresponding mean value for a single stochastic 
system. 

Let us compare the quasi-uniform approximation with the quasi-classical approx- 
imation which is obtained by tending h — > 0. Let us go to the limit h — > in the 
action (4. 15)- (4. 18). One obtains 

AdciU, $]=J [-m^fJJi - eAij 1 - fd^}d 4 x (5.20) 

where 

$ = tup (5.21) 

The term —hfdi(p of Eq. (4.16) can be saved at h — > by means of the substitution 
tp — > $ = Hip. But the last term of C c i cannot be conserved by a like substitution, 
because of the restriction (4.20). For instance, due to Eq.(4.20) the substitution 
£ — > w = h 1 ^ 2 ^ leads to w = at h — > 0. 

The action (5.20) describes only a part of extremals (solutions) of the action 
(5.13), namely that part of them which does not contain spin variables (3.23) and 
describes a potential solution, where the momentum (2.12) forms a potential vector 
field. 

Pi = (5.22) 

The quasi-uniform approximation (5.13) obtained by a proper choice of rather 
smooth initial conditions is more realistic, than the quasi-classical approximation 
(5.20) obtained in the limit H — > 0, because in reality nobody can change the quan- 
tum constant h. Thus, the quasi-uniform approximation is more preferable as a 
classical approximation of the dynamic system Sd- 

Thus, a usage of the quasi-uniform approximation permits to separate out a 
classical part of the action (4.15)-(4.18). The quantum part C q = C q \ + C q 2 of the 
action (4.15) is determined by Eqs.(4.17), (4.18). C q \ contains a specific variable 
k which can be treated neither as a current, nor as a Lagrangian coordinate £. 
Suppressing k (i.e. setting k = 0), the term C q \ vanishes. The term C q2 is a quantum 
term, in general, although it contains a classical part which can be separated out by 
means of introduction of specific quantum variables. 

Let us introduce new variables 

k 1 = q l = —L + Pfl =— f + pr =, i = 0,l,2,3; (5.23) 
VMP + J s fs) V(j s + Pf s )(js + Pfs) 



P = VJJh 

using the Lagrangian multipliers Aj. Then C q2 is substituted by 
C' q2 = hpe lklm (diK l )K t u m + !lfe mm q i j k d s q l is m + 
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+ W + pf - ^V2p(P + J s / s )] (5-24) 
p=y/jijt, u m = [C-f m (Cfs)], m = 0,1,2,3 (5.25) 

where <9^ is defined by Eq.(5.4). Then variation of the action with respect to k % 
leads to the dynamic equations 

5A" 

- \ iy /2p(p + j°f a ) + = 0, (5.26) 

JQj&A = j hpe iklm (d k ± K l )K l u m d A x (5.27) 

Resolving the equation (5.26) with respect to Aj and substituting the A« into 
Eq.(5.24), one obtains instead of C' q2 

C g3 = UpSM^d^ - q l d k ^ + ^d\q l )v m + - h - l \ .. :^J-i)JrC (5.28) 

Z P\P + J Jj) 

Now the action has the form 

A D \j, cp, k, C, «*] = J {Cd + £ g i + C q3 )d 4 x (5.29) 

Dynamic equations generated by actions (3.1), (4.15) and (5.29) are equivalent. 
The quantities £ Q cannot be treated as Lagrangian coordinates (numbers labelling 
systems of the ensemble), because they are not constant along world lines of parti- 
cles, and relations (2.10) do not take place. According to Eq.(3.23) the quantities 
C should be treated as some functions of the spin S = {2S 23 >°, 2S 31 ' , 2S 12 > } = 
{S\S 2 ,S 3 }. 

The variables k, k\ i = 0, 1, 2, 3 are special quantum variables which are respon- 
sible for quantum effects described by the Dirac equation. Indeed, let us set n = 0, 
K l = 0, i = 0, 1,2,3 in the action (5.29). Then the action (5.29) turns to Eq.(5.13) 
with g a {$,) defined by Eq.(4.5). The action (5.29) generates the dynamic equation 



SA D 



he iklm {pu m di(q l - k 1 ) + df[pv m {q l - k 1 ))} = (5.30) 



where the operator d*^ is defined by the relation 

df u = d k u-d s {^u) (5.31) 

Resolving Eq. fl5.30p with respect to k 1 and substituting the k % into Eq.(5.29), one 
returns to the action (4.15)-(4.18). The fact that the solution (5.23) of Eq.(5.30) 
is not sole is of no importance, because Eq.(5.28) reduces to Eq.(4.18) by virtue of 
Eq.(5.30). Indeed, convoluting Eq.(5.30) with k % and using the obtained relation for 
eliminating k 1 from Eq.(5.24), one obtains Eq.(4.18). 
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Let us note that k 1 are not rigorous dynamic variables, because the dynamic 
equations (5.30) for k 1 contain derivatives only along spacelike directions orthogonal 
to f. Rather the introduction of k 1 is an invariant (with respect to a change of 
variables) way of separating out the classical part of the action. 

The concentrated dynamic system Sd c has eleven degrees of freedom. It is asso- 
ciated with the distributed dynamic system Sd- It is described by the action 

A dM - /(-™^-eV + *^| + S ^|^^ (5.32) 

with non-relativistically invariant Lagrangian. It is a very surprising fact which 
needs a special investigation. 



6 Relativistical Invariance 

At our consideration of the relativistic invariance of the Dirac equation written in 
hydrodynamic terms we shall follow the approach of Anderson [6] with the modi- 
fication that the definition of the relativistic covariance is provided by an explicit 
reference to the quantities with respect to which the dynamic equations are relativis- 
tically covariant. Let us consider a simple example which is relevant to the Dirac 
equation. 

One considers a system of differential equations consisting of the Maxwell equa- 
tions for the electromagnetic tensor F tk in some inertial coordinates x 

d k F lk = 47iJ\ e Mm g^d,F kl = (6.1) 

and equations 

™^[(Wry - ^(/^"V W] = eF u g lk (q)q k ; i = 0, 1,2,3 (6.2) 

■k d Q k 

<> s * 

where q % = q l (r), i = 0, 1, 2, 3 describe coordinates of a pointlike charged particle as 
functions of a parameter r, Zj is a constant timelike unit vector, 

9%h = i; (6.3) 

and the speed of the light c — 1. 

This system of equations is relativistically covariant with respect to quantities 
q\ F lk , J\ li, g ik , i.e. it does not change its form at any infinitesimal Lorentz 
transformation 

x* -> x* = x* + w* k x k + o(w); u4 = £*W; u ik = -u ki (6.4) 



19 



which is accompanied by corresponding transformation of quantities q\ F lk , J\ li, 

9ik, 

q\r) ~q\r) = ^? fc (r) = ?* + ^ + o(w) (6.5) 

F lk (x) ->■ F lfc (5) = ^ (x) = + c^F< fe + c^F l/ + o(u) (6.6) 

ji {x ) ^ = -^J k {x) = .r + u\J k + o(u>) (6.7) 

h -)• £ = ^r-/fc = /i + k + o(u), uf = g kl u u (6.8) 

9ik(x) ->■ ^ fe (x) = -Q^i^~k9iA x ) = 9ik + uJi9ik + wigu + o(u) (6.9) 

The reference to the quantities q\ F lk , J 1 , li, means that all these quanti- 
ties are considered as formal dependent variables, when one compares the form of 
dynamic equations written in two different coordinate systems. For instance, if a 
reference to J 1 is omitted in the formulation of the relativistic covariance, it means 
that J 1 are considered as some functions of the coordinates x. If J 1 ^ 0, then 
according to Eq.(6.7) J 1 and J 1 are different functions of the arguments x and x 
respectively, and the first equation (6.1) has different form in different coordinate 
systems. In other words, the dynamic equations (6.1)-(6.2) are not relativistically 
covariant with respect to quantities q\ F lk , U, g ik . Thus, for the relativistic covari- 
ance it is important both the laws of transformation (6.5)-(6.9) and how each of 
quantities is considered as a formal variable, or as some function of coordinates. 

Following Anderson [6] we divide the quantities q\ F tk , J\ U, g ik into two parts: 
dynamical objects (variables) q l , F lk and absolute objects J 1 , U, g^. By definition 
of absolute objects they have the same value for all solutions of the dynamic equa- 
tions, whereas dynamic variables are different, in general, for different solutions. If 
the dynamic equations are written in the relativistically covariant form, their sym- 
metry group (and a compatibility with the relativity principles) is determined by 
the symmetry group of the absolute objects J\ li, git- 
Let for simplicity J 1 = 0. A symmetry group of the constant timelike vector 
li is a group of rotations in the 3-plane orthogonal to the vector U. The Lorentz 
group is a symmetry group of the metric tensor g ik =diag {1,-1,— 1,-1}. Thus, the 
symmetry group of all absolute objects /«, g^, J 1 = is a subgroup of the Lorentz 
group (rotations in the 3-plane orthogonal to li). As far as the symmetry group is a 
subgroup of the Lorentz group and does not coincide with it, the system of equations 
(6.1)-(6.2) is non-relativistic (incompatible with the relativity principles). 

Of course, the compatibility with the relativity principles does not depend on the 
fact with respect to which quantities the relativistic covariance is considered. For 
instance, let us consider a covariance of Eqs. (6.1), (6.2) with respect to quantities 
q l , F lk , J 1 = 0. It means that now k are to be considered as functions of x (in 
the given case these functions are constants), because a reference to U as a formal 
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variables is absent. After the transformations (6.4)-(6.9) the equation (6.2) takes 
the form 

d u ~ dq k x dc£ I 1 r dq k 2 dq s dq s ~, t dq k 

Here U are considered as functions of x. But U are other functions of x, than U of 
x (other numerical constants 4 = Ijdx^ /dx k instead of 4), and equations (6.2) and 

(6.10) have different forms with respect to quantities q l , F lk , J 1 = 0. It means that 
Eq.(6.2) is not relativistically covariant with respect to q l , F th , J 1 = 0, although it 
is relativistically covariant with respect to q l , F tk , Zj, J 1 = 0. 

Setting U = {1, 0, 0, 0}, t = q°(r) in Eq.(6.2), one obtains 

m % =eF o + eF ^ /3 ' < = « = 1,2,3; 
mddq a dq a dq a 

2dt [ dt dt } ~ 6 -°dt ' *" U - (b ' iij 
It is easy to see that this equation describes a non-relativistic motion of a charged 
particle in a given electromagnetic field F lk . The fact that the equations (6.2) or 

(6.11) are non-relativistic is connected with the space-time split into space and time 
that is characteristic for Newtonian mechanics. This space-time split is described 
in different ways in Eqs. (6.2) and (6.11). It is described by the constant timelike 
vector li in Eq.(6.2). In the equation (6.11) the space-time split is described by a 
special choice of the coordinate system whose time axis is directed along the vector 

r. 

Thus, a relativistic covariance in itself does not mean a compatibility with the 
special relativity principles. It is important with respect to which quantities the 
dynamic equations are relativistically covariant 

There is something like that in the case of the Dirac equation. The Dirac equa- 
tion (1.1) is relativistically covariant with respect to variables ip, A iy which are trans- 
formed at the Lorentz transformation (6.4) as follows 

~ 0x k 

Ai(x) -)• Mx) = QpMx) = A i + + o(u) (6.12) 

il>(x) ->■ i)(x) = exp(^i ik u ik )i){x) + o(u), (6.13) 

Y^T = i\ * = 0,1,2,3. 

A direct physical meaning of the variables ip is unclear. Only such quantities as the 
current f and spin (1.6) have a direct physical meaning. The action (4.15)-(4.18) is 
a result of the Dirac action (3.1) transformation to variables j\ S 1 having a physical 
meaning. 

Let us replace the relation (3.23) by the relation 

^ =A ,-i[5'- Z = 0,1,2,3; (6.14) 

P + J S Js 
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where the vector /' is defined by Eq.(4.7). According to relations (1.4) the equation 
(6.14) is equivalent to Eq.(3.23), if the first relation (4.7) takes place. 

Let f k be transformed as a vector, then according to Eq.(6.14) £ l is a pseudo- 
vector, because j l is a vector, and S l is a pseudo- vector. Substituting Eq.(6.14) into 
Eqs.(5.29), (4.16)-(4.17), (5.28), one obtains the action for the Dirac equation in 
terms of variables j, S, k, (p, k 1 

A D \j, S, k, <p, k<] = J (C' d + C ql + C q3 )d 4 x (6.15) 

C' cl = -mp - ej% - hfidtf + e mjkl ^d t p k f m z l ) (6.16) 

where C ql and C qS are defined by Eqs.(4.17), (5.28), p? is defined by Eqs. (4.12), 
(4.11), (6.14), and z/ is defined by Eq.(4.11) 



v 



w k - J ,7 (6.18) 



J .../, f™ m fn 
P + ffs 

For obtaining dynamic equations the action (6.15) is varied with respect to variables 
f, S\ k\ ip, k under additional constraints (1.4). Vector f l and pseudo-vector z % 
are constant and satisfy the conditions (4.8). 

The obtained dynamic equations are equivalent to the Dirac equation (1.1). They 
are relativistically covariant with respect to quantities j l , S l , k 1 , ip, k, f l , z l . 
But they are not relativistically covariant, generally, with respect to variables j l , 
S\ k\ ip, k Ai. It is easy to verify, setting values (4.7) of f\ z % into the action 
(6.15)-(6.18). 

Thus, a relativistical covariance of dynamic equations described in terms of A^, 
j\ S\ k 1 , p, k depends essentially on introduction of additional quantities f\z\ 
These quantities are constant. They arise from the zero divisor (3.5) which can be 
written in the covariant form 

n = l - (i + fa 1 ) \ (i - iiJ lZkl lk ) (6.19) 

where /' is a constant timelike unit vector, and z l is a constant unit pseudo-vector. 

It is easy to verify that under the conditions (4.8) the factors of the matrix 
(6.19) commute, and each factor is a zero divisor. Using designations (3.2), one 
obtains that under conditions (4.7) the expression (6.19) coincides with Eq.(3.5). 
Thus, variables f l ,z l are present in a hidden form inside the zero divisor II of the 
expression (3.4) for the wave function. When one says that the Dirac equation 
(1.1) is relativistically covariant with respect to variables Ai,ip, one implies also the 
absolute objects f l ,z l , hidden inside ip. It is impossible to construct a proper zero 
divisor without using absolute objects. 

Indeed, if all components of the wave function ip are given, the relations (1.2) 
determine the quantities j l and S l uniquely. The quantities <p, k, are determined to 
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within an additive constant and to within an additive (27r)fold constant respectively. 
The quantities do not depend on the form of the zero divisor II. If the quantities 
j l , S l , <f, k are given, the wave function components are determined by the relation 
(3.25) with the zero divisor determined by Eq.(3.5). The wave function components 
depend on parameters determining the form of the zero divisor. In other words, 
the transformation from the variables j l , S l , ip, k to the wave function if) is not 
single- valued. Substitution of from Eq.(3.4) into Eq.(l.l) leads to appearance of 
constant parameters f k , z k in dynamic equations, provided they are written in the 
relativistically covariant form. 

A common continuous symmetry group of both vectors f k and z k is a group 
of rotation inside the 2-plane orthogonal to the bivector f k z l — f l z k . The Lorentz 
group is not a symmetry group of the dynamic equations. If the peculiar directions 
determined by parameters f k , z k are not fictitious, then the dynamic equations 
generated by the action (6.15) are incompatible with the relativity principles. 

The variable f l describes the space-time split into space and time, what is char- 
acteristic for the Newtonian mechanics. In this sense the Dirac equation is not more 
relativistic, than the non-relativistic equation (6.2). 



7 Concentrated Dynamic System Associated with 
the Dirac Equation 

Let us return to an investigation of the dynamic system Sd c - The action (5.32) can 
be presented in the relativistically covariant form as follows 

■ • k 



A dc [x, f] = J {-m^/¥±i - eAiX 1 - 



_|_ ^ ZjklmX X f £ }cIt (7 1) 

2yJx s X s (x l fi + ^/¥x~i) 

where z l and /* satisfy Eqs.(4.8), and the dot means a differentiation with respect 
to invariant parameter r. £ l is a unit spacelike vector orthogonal to the vector /\ 
It can be presented in the form 

e = , t ~ fr]Sfs (7.2) 

V ~V l Vi + (V l fi) 2 

where r\ % is an arbitrary spacelike vector. 
It is easy to verify that 

se = (g + ezM-f l fkW (73) 

V-v s Vs + (v s f s ) 2 
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Then varying Eq.(7.1) with respect to C and using that 8rf is an arbitrary vector, 
one obtains dynamic equations 

^r(*i + KM - ffk) = + tt* - ffk) =0, k = 0, 1, 2, 3 (7.4) 

There are only two independent dynamic equations among the equations (7.4), be- 
cause a convolution of Eq.(7.4) with vectors £ fc and f k converts them into identities. 
The action (7.1) is invariant with respect to a change of the invariant parameter 

T 

t^t = f(r) (7.5) 

where / is an arbitrary monotone function. In particular, r can be chosen in such 
a way that 

x% = 1 (7.6) 

In this case r can be interpreted as a proper time. 

Let us written down the dynamic equations in the case of absence of the elec- 
tromagnetic field Ai = 0. 

5Ad ° ~ d -{^ ~ he tklm [Qfx k C + ±j*±(Qf'e»)]- 



5x i (h ^x s x s ' ' " 2 dr 

— -ZjklmJ X X I; — 



^e jklm fx^ k C^} = 0, i = 0, 1, 2, 3 (7.7) 



se : -{e sklm [pftz m + \z k ^( p f l n + \Q##r\+ 

1 k BP 

2 W'rt + ?ti ~ /Vi) = ( 7 - 8 ) 
where P and Q are defined by the relations 

Q = Q(x) = [v^(i7i + V^)]^ (7-9) 

P = P(0 = (1 - z 1 ^)- 1 (7.10) 
Equation (7.7) is integrated in the form 

' He Mm [Qf l S*C + \i k ^-{Qf l n\ ~ ^ jklm f l x^x k r^ = p t (7.11) 



y/¥¥ a ■ — ^ J ~ ^ • 2 dr^ J * n 2 * dx* 

Pi = const, i — 0,1, 2, 3; 

Let r be the proper time, and the condition (7.6) takes place. Convoluting 
Eq.(7.11) with /' and eliminating the last term in lhs of Eq.(7.11), one obtains 

- e iklm [x k fr + \x k fT ~ ^lf s J^ = - fif S )( u s - is)^-^, (7.12) 

A = — , iii = —, % — 0,1,2,3 
m m 
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One can see from Eqs.(7.8), (7.12) that only directions orthogonal to the timelike 
vector /* are essential. 

Let us choose /' in the form (4.7). Then using designations 

x* = {x°,x}, f = {0,£}, ^ = {0,z}, ^ = { M °,u}, (7.13) 

Eqs. (7.8), (7.12) reduce to the form 

K + , (Z *> g] x z + ^ XZ K + ii±4-x xx = C« (7.14) 

L S 2(l + z£) SJ 2(1 + z£) 2(1 + a: ) ^ V ; 

d t x \ xx£ y'l + x Q . . . ,., r > 

"tt( 75= ^ x € + 9 7== = \ x - u (7.15) 
ctr ^l + x 2y/l + x A 

where C is some indefinite function of r. 
According to Eqs. (7.2), (4.7) 

£ 2 = l, « = 0. (7.16) 
Using these relations, the equation (7. 14) reduces to the form (see Appendix B) 

• (x x x) 



l + x 



x £ (7.17) 



which does not contain the vector z. It means that z determines a fictitious direction 
in the space-time. Note that z in the action (4.15) for the system S& is fictitious 
also, because the term containing z is the same in the actions (4.15) and (5.13) for 
Sd and So q u respectively. 

Eq.(7.17) can be written in the relativistically covariant form 

= TTfJ* (7 - 18) 

The system of Eqs. (7. 12), (7.18) is relativistically covariant with respect to the 
quantities x\ £\ p iy f % considered as 4-vectors (x\ p iy f l ) and pseudo-vectors (£*). 
But it is not relativistically covariant with respect to the dynamic variables x l , £ l ,Pi- 
Under the Lorentz transformation (6.4) the set Sf of all solutions {x l ,l;\pi} of 
Eqs. (7. 12), (7.18) at fixed parameters f l is transformed into another set Sf of so- 
lutions {x\^\pi}. But the set S = {Sf} of all sets Sf is transformed into itself, 
provided /' satisfies Eq.(4.8). 

Formally the integration constants pi can be considered as some parameters of 
the system of differential equations (7.18), (7.12). At the same time the parameters 
Pi, f l of Eqs. (7.18), (7.12) can be considered as some integration constants of some 
system V of differential equations which is relativistically covariant with respect to 
its dynamic variables. The set of all solutions of V can be obtained as a result of 
the Lorentz transformation (6.4) of the set S p of all solutions with fixed values of 
Pi, or as that of the set Sf of all solutions with fixed values of f l . 
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It is sufficient to investigate solutions Sf of equations (7.12), (7.18) at fL, = 
{1, 0, 0, 0}, because solutions at other values of f 1 can be obtained from solutions at 
f l = fL, by means of a proper Lorentz transformation. 

Eq.(7.17) describes a rotation of the unit 3- vector £ around the vector x x x 
with the angular frequency 

= x^x 

5 l + x 1 ; 

Eliminating £ by means of a substitution of Eq.(7.17) into Eq.(7.15), and introducing 
a new variable y instead of x 



y =_L=, x = yv /2 + y 2 , x° = VlT^, (7.20) 



Eq.(7.15) reduces to the form 



Ayx(£ + —^=y) = (y v / 2 + y^-u) v / 2T^ (7.21) 

2 v 2 + y 

Let us introduce a dimensionless variable w which is supposed to be small 

w = y-b, b = u(l + Vl + u 2 r 1/2 , | w |<1 (7.22) 
Then Eq.(7.21) takes the form 

Aw x a= |Gw + b(wb)] + 0(w 2 ), (7.23) 

It follows from Eq.(7.23) that a is orthogonal to rhs of Eq.(7.23), or 

wd = 0, d =Ga + (ba)b (7.25) 

Characteristic frequency for w is of the order of A -1 , whereas that for £ is of the 
order of A -1 | w |<C A -1 . It means that, integrating Eq.(7.23), one can consider £ 
as a constant vector. 

If £ || b, w is orthogonal to b, and Eq. (7.23) reduces to the form 

G 

Aw = — a x w + 0(w 2 ), w 2 = const (7.26) 

which describes a rotation of the vector w with the angular velocity Ga/a 2 . 

If £ x b t^O, let us choose orthonormal vectors e 1? e 2 , e 3 in such a way that ei, e 2 
lie on the plane V spanned by vectors £ and b. Let the vector e x be chosen on V in 
such a way that ei || d. Then one has the following decomposition of vectors w, a, b 

w = w 2 e 2 + w 3 e 3 , a = aiei + a 2 e 2 , h = + b 2 e 2 (7.27) 
According to Eqs.(7.23), (7.25), (7.27) the first equation (7.25) reduces to the form 
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w 2 [ai&i&2 + a 2 {G + b 2 2 )] = 0. (7.28) 
Substituting Eqs. (7.27) into Eq.(7.26), one obtains 

-Xw 3 a 2 = b 1 b 2 w 2 + 0(w 2 ), 

\w 3ai = (G + b 2 2 )w 2 + 0(w 2 ), (7.29) 

— \w 2 ai = G1V3 + 0(w 2 ), 

Due to Eq.(7.28) two first equations (7.29) are equivalent. The system (7.29) can 
be easily solved. The solution has the form. 

w = A[e 2 cos + (1 + b 2 2 /G) 1/2 e 3 cos 0] + 0(A 2 ) (7.30) 

where 

ei=d/|d|, e 2 = [1 - (ei^) 2 ]" 1/2 ei x £, e 3 = e x x e 2 , 

= fi x (r-r ( o)) + (O ) (7.31) 

n x = ^-[l + bj/G] 1 / 2 , ai = (a ei ), b 2 = (be 2 ) (7.32) 
Using Eqs. (7.20), (7.22), one obtains 

x = u + A[(Ge 2 + ^6) cos0 + (G + 6 2 ) 1/2 e 3 sin0] + 0(A 2 ) 

x° = VT+i 2 (7.33) 

x = x {0) + u(t - r (0) ) + Afi; 1 [(Ge 2 + ^6) sin - (G + 6 2 ) 1/2 e 3 cos 0] + 0(A 2 ) 

Gr 



T 

J 0) + | aAW)^' (7.34) 

T (0) 

A, 0(0-), x° ^, X( ), u, T( ) are independent integration constants, a, b, G and Q x are 
determined by relations (7.22), (7.24), (7.32) through these constants. 

According to Eq.(7.19) the angular velocity Q$ of the vector £ rotation reduces 
to the form 



^ = — — — t^j- [ — u x (Ge 2 + ^6) sin0 - (1 + 6 2 G) 1/2 u x e 3 cos0] + 0(A 2 ) (7.35) 

oscillates rapidly around the mean value < J\ >~ 0(a 2 ). It means that 

Z = Z (0) +O(A 2 ) (7.36) 

£(0) = const > £(0))£(0) = 1 ( 7 - 37 ) 
27 



In order to show that f l is not a fictitious parameter, let us consider the exact 
solution of equations (7.12), (7.18) at /* defined by Eq.(4.7) . 



{Vl + a 2 r, -aft -1 sin 0, aft -1 cos 0, 0}; f = {0,0,0, 1} 



1 + Vl+V 



a 2 

u l = {1,0,0,0}, n = - ' v ^ ' , = fir + (o) , a = const (7.38) 

If f 1 is a fictitious, the system of equations (7.12), (7.18) is relativistically covari- 

ant with respect to the vectors x\£\ u \ an d the vectors x\ | , u\ obtained from 
Eqs.(7.38) by means of a Lorentz transformation are to form a solution of Eqs.(7.12), 
(7-18). 

In the coordinate system moving with the velocity V = tanh-$ in the direction 
of the axis x 3 the vectors (7.38) have the form 

x 1 = {tV 1 + a 2 cosh — af2~ 1 sin0, af2~ 1 cos0, ry/l + a 2 sinh-^}; 

|* = {sinhtf, 0, 0, coshtf}; u* = {coshtf, 0, 0, sinh^}; (7.39) 



1 + y/TT^ 2 " 

sZ = = s 2r + 0(0-) , a = const 

It is easy to verify that (7.39) is not a solution of Eqs.(7.12), (7.18), if $ ^ and 
/' is determined by the equation (4.7). (Of course, Eq.(7.39) is a solution of Eqs. 
(7.12), (7.18), if /* = /* = {cosh^,0,0,sinh^}). It means that the vector /* is an 
essential parameter, and the system of equations (7.12), (7.18) is incompatible with 
the special relativity principle. 

The fact that z l is fictitious, but f l is not, is explained, apparently, by incomplete 
symmetry of the Dirac equation with respect to the time and the space. The matrix 
7° is used for constructing the Dirac conjugate spinor ip = -0*7° • It separates out 
the matrix 7 among matrices 7* 

A space-time split generated by 7 appears in the Space-Time Algebra (STA) 
suggested by Hestenes [7]. STA is a kind of Clifford Algebra describing space-time 
properties [8]. This space-time split is connected with a use of the matrix 7 . 

In the case, when the electromagnetic field does not vanish, the quantities Pi 
determined by Eq.(7.11) are not constant. They satisfy the equations 

Pi = mui = eF ik (x)x k , F ik = diA k - d k A, L . (7.40) 

There are only three independent equations among them 

mix = eEx° + e[x x H], E = E(x°,x), H = H(x°,x) (7.41) 

where E = {E a } = {-F a0 }, a = 1,2,3; H = {H a } = {-\e a ^F^}, a = 1,2,3. 
The system of equations (7.6), (7.16), (7.17), (7.20), (7.41) is a complete system of 
dynamic equations for the dynamic variables £ ,x, x°, u. 

For not too strong electromagnetic field, when the Larmor frequency is much less, 
than il x (eH/m <C £l x ~ m/h, i.e. if H <g; 10 12 G), u in Eq.(7.21) can be considered 
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approximately as a constant. Then the world line x l = x l (r), i = 0, 1, 2, 3 describing 
a solution of the system (7.6), (7.16), (7.17), (7.20), (7.21), (7.41) spirals round the 
world line X % = X 1 {t) which describes a motion of the guiding center and satisfies 
the equation 

mX l = eF l fc (X)X fc , i = 0, 1, 2, 3 (7.42) 

In general, a charged particle moving in such a way must intensively radiate 
electromagnetic waves. As a result its world line is to approach to the world line of 
the guiding center. 

Such a helical motion of the electron associates with the classical model [9] of 
the Schrodinger zitterbewegung [10]. Another approaches to the interpretation of 
the zitterbewegung can be found in papers [11-14] and references therein. 

8 Discussion 

It seems that pi and x l should be interpreted respectively as a 4-momentum and 
4-velocity of the system Sd c - But then the system Sd c can be hardly interpreted as 
a pointlike classical particle. There are two reasons. First, pi and x l of a pointlike 
particle are connected by an algebraic relation of the kind of Pi = mx^. For the 
dynamic system Sdc a distinction between y % = x l (\ + x°) ^ 2 and u l = p l /m is 
described by the dynamic variable w l which is introduced by Eq.(7.22). According 
to Eq.(7.21) this distinction is proportional to the quantum constant h. Second, it is 
very difficult to understand, why the world line of a free particle is a helix spiralling 
round the straight line 

tf = £(T-T (0 )) + Xfo. (8.1) 

even in an absence of the electromagnetic field. 

All this can be explained by means of the supposition that x l describes an observ- 
able part of a complicated bound system whose center of inertia X % moves according 
to Eq.(8.1) in the case of Ai = and according to Eq.(7.42) in the case of not too 
strong electromagnetic field. But under such an interpretation the system Sdc is 
not a pointlike particle with some inner degrees of freedom, but a dynamic system 
consisting of a few particles interacting through a distance. But such an interaction 
through a distance is incompatible with the relativity principle. 

One can try to save the compatibility with the relativity by a consideration of 
some distributed classical system S^ instead of the concentrated dynamic system 
Sdc- The system Sd c arises as an element of the statistical ensemble S^qu- But an 
one-parameter statistical ensemble of dynamic systems Sd c can be also considered 
as an element of the statistical ensemble So q u- For instance, let us consider a set 
Sc/, of systems Sd c , having all similar integration constants {a, x^, X( ), u, T( ), £(o)} 
except of (f), y Let values of <f>, \ be distributed uniformly among the systems Sd c - 
(For instance, in Eq.(7.38) x l is considered as a function of two variables r and 
0( O ) with fixed parameter a). Then S^ can be considered as a ring. World lines 
of the ring particles form a two-dimensional world tube in the space-time. Such 
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a construction is more symmetric, than helically moving particle. Besides such a 
ring does not radiate electromagnetic waves. The statistical ensemble Sdqu can be 
considered as consisted of rings S^ which are classical distributed dynamic systems. 
Unfortunately, rings S^ and their two-dimensional world tubes are incompatible 
with the concept of a pointlike particle. 

One can try to associate the world tubes with the space-time properties (but 
not with the structure of the dynamic system S^). In this case it is necessary to 
admit that the distortion of the space-time does not vanish, i.e. the characteristic 
geometric structures of the space-time are three-dimensional world-tubes (but not 
straight lines) [2,15]. Then one can hope to remove incompatibility between the 
concept of a pointlike particle and the cylindrical space-time structures generated 
by the Dirac equation. Note that in the distorted space-time the momentum of 
a particle and its velocity are independent quantities which are connected by no 
algebraic relation. 

Thus, the well known dynamic system Sd described by the Dirac equation has 
been investigated simply as a dynamic system without using any additional suppo- 
sitions. The obtained results can be summed as follows. 

1. Describing the dynamic system Sd in terms of hydrodynamic variables, one 
discovers that in the quasi-uniform approximation the Sd turns to a statistical 
ensemble Soqu of classical systems Sdc- It means that, in general, the dynamic 
system Sd is a statistical ensemble of stochastic systems, and it permits to separate 
out a classical part C,Dqu of the total Lagrangian Cd- 

2. It is possible to separate out the classical part C-Dqu of the Lagrangian by 
means of special quantum variables k, k % responsible for quantum effects. Suppress- 
ing the quantum variables (setting k, k 1 = 0), one transforms the Lagrangian Cd of 
the system Sd into its classical part Cd qu - 

3. The transformation of the Dirac equation to hydrodynamic variables separates 
out a peculiar direction in the space-time. This direction is described by the timelike 
unit vector f k . Existence of a preferred direction in the space-time is incompatible 
with the special relativity principle. It means that the Dirac equation is incompatible 
with the special relativity principle and needs a modification for such a compatibility. 

4. The dynamic system Sd c is a classical analog of the quantum Dirac electron. 
Sdc cannot be interpreted as a pointlike charged particle in the Minkowski space- 
time. Rather it is either a complicated non-local dynamic system in the Minkowski 
space-time, or a pointlike particle, but in a non-Riemannian (distorted) space-time. 



Classical Ensemble Described in Lagrangian and Eulerian Coordinates 

Let us show equivalency of two form (2.1) and (2.2) of the action for the statistical 
ensemble of classical systems. Let us introduce the time Lagrangian coordinate £ , 
supposing that t = t(£ ), x = {t, x} = £ = {£o>£}- Then the action (2.1) can 
be rewritten in the form 



APPENDIX A 




(Al) 
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C(x, x) = L(x°, x, ±/x°)x°. (A2) 

According to Eq.(A.2) the Lagrangian C(x,x) is a first order homogeneous function 
of x, i.e. 

C(x, ax) — aC(x, x), (A3) 

where a is an arbitrary parameter. 

Let the variational problem with the functional (A.l) take place. Let us show 
that it is equivalent to the variational problem with the functional (2.2). 

Considering £ = i = 0,1,... n as functions of x — {x 1 }, % = 0,1,... n, 

the variational problem with the action (A.l) can be formulated as a variational 
problem with the action 

A[£] = [ C(x, ^-)d n+l x, £ = (A4) 



where J is a Jacobian 



and 



t _ ^(^Qi ^1) • • • Cn) _ j . II p || 

ii, k = dijdx k , i,k = 0,l,...n (A5) 

^■■y ^J, < = 0,l,...n (A6) 

Deriving Eq.(A.4) from Eq.(A.l), the relations (A. 3), (A. 6) were used. 

In the variational problem with the action (A. 4) the Lagrangian coordinates 
£ = « = 0, 1, . . .ra are dynamic variables considered as functions of coordinates 
x = {x 1 }, i — 0, 1, ... n. 

Let us use designations 

(97 

/ = -— , i = 0,l,...n (A7) 

Adding constraints (A. 7) to the action (A. 4), one does not change the variational 
problem, because Eq.(A.4) does not contain variables j\ and Eq. (A. 7) is not a 
constraint in reality. Let us introduce Eq.(A.7) into the action (A. 4) by means 
of Lagrange multpliers p = {pi}, % = 0,1,... n. j l and Pi, i = 0,1,... n mean 
respectively a current density and a momentum density in the space {£, x}. 
Then one obtains 

/B T 
[L(xJ)-iHU i --M-)]dx n+1 . (A8) 

Dynamic equations generated by the action (A. 8) have the form 
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£ = -* fe a& ) = 0, ° = 1 ' 2 -"" i (A9) 

f = §-B = 0, i = 0,l,...« (-4.10) 

Variation with respect to leads to Eq.(A.7). 
Using identities 

d k { d Q ) = 0, i,Z = 0,1,... n, d k j k = d k -j— = 0, (All) 

£ =*f^_3™ i,fc,Z = 0,l,...n, (A12) 



p fc = — = ^ + ^ii)dk^p\ k = 0, 1, . . .n, (A13) 



and eliminating the variables £, one can reduce equations (A. 7), (A. 9), (A. 10) to 
the conventional hydrodynamic form (2.9), (2.14) 

Instead of elimination of variables £ Eqs. (A. 9), (A. 10) can be integrated in the 
form 

dC{x, j 
dj 

where g 13 ^), (3 = 1,2, ...n are arbitrary functions of the argument £ = {£ a }, 
a — 1, 2, . . .n. y9 is some new variable. /I is the Planck constant which is introduced, 
for variables if, £ = « = 1, 2, . . .n and the functions g' 3 were dimensionless. 

In this context h has no quantum meaning. 

Indeed, let us substitute relation (A. 13) into Eq.(A.9) and use the first Eq.(A.ll), 
antisymmetry of d 2 J/d£ 0l d£ ak and symmetry of d k diip with respect to indices I, k. 
Then one obtains 

dk[ g'WZu] = 0, a = 1, 2, ... n. (A14) 

Now using the identity 

d 2 J . _ &7 



^ ^ ^,i = -^^» a,0 = l,2...n, fc = 0,l,...n, (A15) 



and the second relation (A. 11), one can reduce Eq.(A.14) to the form 



Due to the identity 



= 0> a = 1 > 2 >-..n. (Ai6) 



(9 7 

= / = 0,l,...n (A17) 



the equation (A. 16) is satisfied for any functions g 13 ^), (3 = 1,2, ...n. Thus, the 
relation (A. 13) is a general solution of the system of equations (A. 9). 
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Let us substitute Eq.(A.13) into the action (A. 8) and use that due to Eq.(A.17) 

dJ -v t = ^rH—y . (-4.18) 



d$ /" d(x°,x\...x«) ■ 

Then the action (A. 8) turns to Eq.(2.2). The term (A. 18) is omitted, because it 
does not give any contribution into dynamic equations. Thus, Eq.(2.2) is a corollary 
of Eq.(2.1). 

Now let the variational problem with the functional (2.2) take place. Corre- 
sponding dynamic equations have the form (2.9)-(2.11). The relations (A. 7) satisfy 
Eqs. (2.9), (2.11) identically, and an usage of them as additional constraints does 
not change the variational problem. Substituting (A. 7) into Eq.(2.2) and using iden- 
tities (A. 17), (A. 18), one obtains Eq.(A.4) which is equivalent to Eq.(2.1). Thus, 
the variational problems with actions (2.1) and (2.2) are equivalent. 

APPENDIX B 

Let us show here that under conditions 

e = i, z 2 = i (B.i 

the equation 

where b is a given vector, and C is an indefinite function of r, reduces to the form 

i = b x £ (B.3 
Let z£ = cos a ^ ±1. Let us introduce an orthonormal basis {e a }, a — 1,2,3 

e 3 = £, e 2 = -, £ = cosa IB A 

sin a 

£ x (£ x z) z - £cosa . „ _ 

ei = e 2 x e 3 = P = * (B.5 

sin a sin a 

Then 

z = ei sin a + e 3 cos a, £ = £ a e a , b = b a e a (B.6 
It follows from (B.I) and (B.6) that 

U = 0, 4 = 0, z 2 = 0. (B.7 

Substituting Eqs. (B.6), (B.7) into Eq.(B.2) and equating coefficients before the basis 
vectors e a , one obtains 

e x : t 2 = -h (B.8 
e 2 : ti = h (B.9 
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„ b 3 2 + cos a ■ 

e 3 : C = -(l + COBa )-^—^^ 2 . (5.10) 

Eq.(B.3) follows from Eqs.(B.7)-(B.9). Eqs.(B.8)-(B.10) determine the indefinite 
function C through the given vector b and the scalar product (z£). 
If z = £, the equation (B.2) reduces to the form 

-£x£-h = C£. (5.11) 

Forming the vector product of Eq.(B.ll) with £ and using (B.7), one obtains 
Eq.(B.3). 

If z = — £, the equation (B.2) becomes indefinite, and it is necessary to use the 
limit £ — z. It leads to the result (B.3). 
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